Profiles of near-resonant population-imbalanced trapped Fermi gases 
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We investigate the density profiles of a partially polarized trapped Fermi gas in the BCS-BEC 
crossover region using mean field theory within the local density approximation. Within this 
approximation the gas is phase separated into concentric shells. We describe how the structure 
of these shells depends upon the polarization and the interaction strength. A comparison with 
experiments yields insight into the possibility of a polarized superfiuid phase. 



INTRODUCTION 

Advances in the trapping and manipulating of degen- 
erate Fermi atoms are attracting intense interest from 
physicists in the fields of condensed matter physics, 
atomic molecular and optical physics, nuclear physics, 
astrophysics, and particle physics. Current and future 
experiments aim to use this highly controlled environ- 
ment to explore many-body phenomena with impact on 
widely varying areas of physics. We discuss the theory 
of one such set of phenomena; the properties of trapped 
partially polarized fcrmionic atoms, where the two spin 
states have different Fermi surfaces. 

The study of Fermi systems with mismatched Fermi 
surfaces began with attempts by Fulde and Ferrell and 
Larkin and Ovchinnikov (FFLO)|lJ to understand mag- 
netized superconductors. More recent work has come 
from researchers in the nuclear physics community who 
are studying superconductivity in nuclear matter and 
quark matter, with a possible application to neutron 
stars or heavy ion collisions |2j. Such calculations have 
taken on new relevance with the possibility of cold gas 
experiments where alkali atoms are trapped in a num- 
ber of distinguishable hyperfine states, with negligible 
spin relaxation. Thus one can produce a two-component 
Fermi gas with arbitrary population imbalance. Using 
magnetic- field driven Feshbach resonances |3j, the inter- 
actions between these atoms can be made large enough 
to drive the system superfiuid. 

Very recently, there have been two experimental stud- 
ies of trapped 6 Li Fermi atoms with mismatched Fermi 
surfaces |4j5| . By analyzing time-of- flight images, Zwier- 
lein et al. jj| found evidence for phase separation between 
regions of equal and unequal population density. Further- 
more, by studying vortices, they were able to monitor 
the evolution of superfluidity as a function of population 
imbalance: finding not only that polarization inhibits 
superfluidity, but that the superfiuid region appears to 
coincide with the region of equal density. The equally 
exciting work of Partridge et al. directly shows phase 
separation through high resolution in-situ images of the 
atomic clouds. 

In this paper, we study the density profile in the entire 
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FIG. 1: Mean- field phase diagram of homogeneous two- 
component Fermi gas in the (a), (b) BEC and (c) BCS 
regimes. Phases only depend on the dimensionless ratio of 
the chemical potentials fi — (/tf + /tj,)/2 and h — (/if — /t|)/2, 
and the energy scale of the interactions, e = h 2 /2ma 2 . Panels 
(a) and (b) differ only by the scale used. 



BCS-BEC crossover region within mean field theory. We 
use a local density approximation (LDA) and compare 
our results with experiments. Our study shows that the 
population imbalanced trapped Fermi gas is generically 
phase separated into concentric shells. Within our ap- 
proximations, each region of space can be in one of several 
phases: unpolarized superfiuid (S), polarized superfiuid 
(PS), normal mixture (M), or fully polarized normal (P). 
The unpolarized superfiuid coincides with the standard 
equal-population superfiuid predicted by the BCS-BEC 
crossover theory |6|, LI la l2l llfll • The polarized superfiuid, 
which in mean-field theory is only found on the BEC 
side of resonance, consists of an interpenetrating gas of 
bosonic molecules and a fully polarized gas of fcrmions 
[ill . The normal mixture and fully polarized normal 
phases both lack superfluidity and are distinguished by 
the presence or absence of the minority species of fermion. 
Due to the extremely small portion of the phase diagram 
where it is expected to appear, we do not consider the 
possibility of an inhomogeneous FFLO phase. As seen 
in a number of recent papers 0, 0, Q] , such a phase 
would appear on the BCS side of resonance as a sub- 
tle structure in the domain wall separating the S and M 
regions. 

We find three different parameter regimes, distin- 
guished by the structure of these shells: (1) BCS regime 
where (fc/a) -1 < 0, (2) crossover-BEC regime where 
< (fc/a) ~ 1) an d (3) deep-BEC regime where 
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(kja)^ 1 > 1. The Fermi wave vector is kf and the scat- 
tering length is a. Depending upon density and polariza- 
tion, regime (1) contains two scenarios - from center to 
edge one finds S-M-P or M-P. In regime (2) one finds S-P 
at small polarizations and S-M-P or M-P for larger. In 
regime (3) one finds S-PS-P for small polarizations and 
PS-P for larger polarizations ^5|. In an appropriately 
defined thermodynamic limit there is a quantum phase 
transition between each of these possibilities as one varies 
the parameters of the system. This behavior should be 
contrasted with the smooth crossover physics found in 
the absence of a population imbalance. 

Our results are consistent with the experiments re- 
ported in Ref. 0], however, due to the expansion pro- 
cedure used in those experiments no quantitative com- 
parison can be made. We find partial agreement with 
the experiments reported in Ref. In particular, for 
sufficiently large polarizations we reproduce the values of 
the axial radius of the superfluid core and the outer edge 
of the gas cloud found in the experiments 0- Our total 
density profiles also closely agree with the experiments. 
However, we find that the spatial structure of the dif- 
ference between up and down spins differs significantly 
from those found in the experiment []|. In particular, we 
show that despite the large ratio of the trap size to the 
coherence length, the experimental data are inconsistent 
with any LDA that assumes a harmonic trap, regardless 
of the equation of state used. 

The most significant open theoretical question at the 
moment is the possibility of a polarized superfluid re- 
gion at unitarity (UPS). As seen from the e — ► behav- 
ior of Fig. db), the mean-field calculation predicts that 
such a region does not exist. Monte Carlo calculations 
suggest that such a region may exist, but are currently 
not conclusive Based upon a comparison with our 
mean-field calculations we argue that the experimental 
measure of phase separation (from analyzing the density 
profiles) is slightly ambiguous and one may be able to 
explain the experiments without recourse to a UPS. On 
the other hand, the poor agreement in the radii at small 
polarizations may support the notion of a UPS. 

Concurrent with the preparation of this manuscript 
several authors [l^L Il7j | presented complementary theo- 
retical studies of the effect of trap potentials on a par- 
tially polarized gas. With the exception of Ref. 11], which 
discusses some qualitative feature of the trapped gas, 
prior theoretical work on superfluidity with mismatched 
Fermi surfaces has been restricted to either homogeneous 
systems or trapped systems in the 

weakly coupling limit |l4| . 

THEORY 

We restrict ourselves to the wide resonance of 6 Li 
where both of the available experiments 0, 0| have been 
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FIG. 2: On-axis mean-field density profiles n(z, p = 0) 
of zero temperature harmonically trapped partially polar- 
ized Fermi gas, illustrating all of the cases described in the 
text. Solid lines show total density n-\ + nj., while dashed 
lines show density differences rij — rij. Left (right) fig- 
ures represent the strong coupling BCS (BEC) regime with 
kfd — —2, (0.5). The Thomas-Fermi radius is defined as 
Rtf = y/(2ef)/(mu$), where e/ = h 2 k 2 f /2m = huj(6N) 1/3 
with average trap frequency ZJ — (tj^u z ) 1//3 . Inset shows col- 
umn density rv c ' = J dpn(z,p) measured in units of [10 9 
cm -2 ]. All graphs use ATf + N j, = 2 x 10 5 atoms in an 
axial symmetric trap with uo z = 2n x 7.2Hz, and u>± = 
2ir x 350-ffz. BCS figures (top to bottom) have polarizations 
P = (ATf — JVi)/(JV T +N k ) = 0.30, 0.74, 0.80 and BEC figures 
have P = 0.33, 0.90. 



carried out, and therefore we do not need to explic- 
itly consider the closed channel of the Feshbach res- 
onance. The fermions of different hyperfine states J 
and I interact through a short-range effective poten- 
tial —u5(ff — r). The system of N = iVj + iVi atoms 
is then described by the Hamiltonian 0, 0, Ej H = 
Hi + H 2 , with Hi = E CT /^t(r)[-^ - Po* + 
U(r)]i/v(r) containing kinetic and trapping energies and 
H-2 = —u J d 3 f^l(f^ip](f^'ipi(r)'ip'i(r) containing interac- 
tions. The field operators ^ a {r) obey the usual fcrmionic 
anticommutation rules, and describe the annihilation of 
a fermion at position r in the hyperfine state a. Pa- 
rameters rn, /ioo- and V(r) — ^m(uj\p 2 + w^z 2 ) are the 
mass, chemical potential, and trapping potential of the 
atomic species a. We introduce a local chemical poten- 
tial fia(r) = /loer — V(r) and treat the system as locally 
homogeneous. Without loss of generality we take f to 
be the majority species of atoms and describe the system 
in terms of the spatially independent chemical potential 
difference h = (p^ — pi)/2 > 0, and the spatially de- 
pendent average chemical potential p(r) = (p^ + pi)/2. 
Using the usual BCS mean-field decoupling, the BCS- 
Bogoliubov excitation spectrum of each species is given 
by E ka {r) = tah+^/(e k - pf + A 2 . Here e fe = h 2 k 2 /2m 
is the kinetic energy, A(r) = u(if;i(r)t(jy(r)) is the local 
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superfluid order parameter, and £j = —1 and ^ = +1. 

At zero temperature, the gap A(r) and the the number 
densities n(f) = n^(f) + n^r) = nj (r) — n±(r) 

satisfy 
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We define ^(r) = 
(±V/i 2 - A 2 +M) 1/2 - 



(£ feT + £ fci )/2, and fc±(f) = 
The ultraviolet divergence asso- 



ciated with the delta- function interaction has been elimi- 
nated 0,00,113 by introducing the effective scattering 
length through -m / {A-nfi 2 a) = v,- 1 ~ J °° d 3 k(2ir)- 3 /2e k . 
Notice that ridir) is nonzero as long as h > A and 
H > -Vh 2 - A 2 . 

At each point in space one finds A by solving Eq. 
at fixed p and h. As a nonlinear equation, there are 
multiple solutions - some of which correspond to energy 
minima, some to saddle points. For example, the Sarma 
state Pii I21J Hl| appears as a saddle point. We take the 
lowest energy solution, producing the phase diagram il- 
lustrated in Fig. n Within the cloud h is uniform, while 
p varies monotonically from the center to edge of the 
cloud. The configurations of local phases, as described in 
the Introduction, are found by following vertical lines in 
the figure. One determines po and h by imposing a con- 
straint on the total number of particles N = J d 3 fn(r) 
and the polarization P = (Aj — N±)/N. Typical density 
profiles are shown in Fig. [21 Evolution of the radii of 
phase boundaries are shown in Fig. |3| 



EXPERIMENTS 

There are profound differences in the full three- 
dimensional density profile of the experimental cloud of 
Ref. and our predictions. These differences are best 
seen by looking at the axial density profile n j>°^ (z) — 
2ir J dpprid(z,p), found by integrating the densities over 
both transverse directions. As shown in Fig. we find 
a monotonic axial profile, while at polarizations above 
P « 0.1, Partridge et al. experimentally see a non mono- 
tonic density difference, with a dip in the center and 
horns on the edges. Despite these horns, our calcula- 
tion of the total axial density n^ a \z) — 2tt J dppn(z,p) 
agrees extremely well with the experimental data (Fig. 
EJ. 

Assuming LDA in a harmonic trap, the density is 
only a function of p(r) = po — m,Lo\p 2 /2 — mcu 2 z 2 /2. 
[As before, p = (mot + Mo|)/2-] One can therefore 
write n { d a) (z) = ^rf(po - muj 2 z z 2 /2), where f(p) = 
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FIG. 3: TOP: Radii of the minority and majority compo- 
nents using experimental parameters from Ref. [fj. Dots and 
crosses are the experimental data for the majority and mi- 
nority components respectively. Lines show our theoretical 
predictions with fc/a = 20.6. The inset shows the total num- 
ber of atoms used for calculation at each P. The two spin 
state radii R±) are separately scaled by Rtfi and Rtfi, 
where Rtf<j is the Thomas- Fermi radius of a noninteracting 
one-component cloud with N a atoms. BOTTOM: Radii of 
phase boundaries in the BCS regime with kfa = — 2. Solid, 
short-dashed, and long-dashed lines are the boundaries of po- 
larized normal, normal mixture, and unpolarized superfluid, 
respectively. Input parameters are the same as those of figure 
|5J as is the definition of Rtf- 
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FIG. 4: Doubly integrated axial density difference ri^\z) = 
2n J dp p[ni(z, p) — ni (z,p)] of zero temperature harmon- 
ically trapped partially polarized Fermi gas in units of 
[10 6 cm _ ]. Figures (a), (b), (c), and (d) represent polariza- 
tion P = 0.01, 0.09, 0.14, and 0.53 respectively. All graphs 
use the experimental parameters from Ref. |f| and the gray 
points are the experimental data. 



f_ dpnd(p)- Since rid > one has f(p) is monotonic 

and n d a \z) must decrease monotonically as z increases 
from 0. In other words the horns seen in the experimental 
density differences are not consistent with the LDA. This 
result is based solely on the local density approximation 
and harmonic trapping - it does not require mean-field 
theory. We caution that this theorem only applies to the 
doubly integrated axial density, and not to the column 
density. 

Even if mean-field theory breaks down, these exper- 
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FIG. 5: Comparison of the axial density n^ a \z) = 
2-rr J dp p{n^{z, p) + n-i(z,p)] in units of [10 6 cm -1 ] with ex- 
perimental data of Ref. [f| • The labels are the same as those 
of Fig. H 



iments should be well described by a local density ap- 
proximation. The only relevant microscopic scale near 
unitary is the Fermi temperature Tf « 400^^, which is 
over 20 times larger than the quantization scale of the 
harmonic trap Tiu>_\_/kB ~ 17nK, which is a characteris- 
tic scale for density variations. We are currently investi- 
gating the possibility that surface tension in the domain 
wall between the polarized and superfluid regions may 
be distorting the shape of the boundary, leading to the 
observed density profiles. 

Further experimental and theoretical work is needed 
to settle this issue. It is particularly important because 
Partridge et al. interpret the appearance of horns as a 
transition between a unitary polarized superfluid and a 
phase separation between superfluid and normal regions. 
Since the LDA predicts that the horns do not exist in a 
phase-separated cloud, we caution strongly against tak- 
ing their disappearance as evidence for a polarized super- 
fluid. 
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